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Key Objectives 


• To study the 29-30 May 2012 deep convective storm observed 
during the Deep Convective Clouds and Chemistry (DC3) 
experiment over Oklahoma, including its: 

- Convective transport of trace gases 

- Associated lightning occurrence and nitrogen oxide (NO x ) production 

• Simulate the observed storm using WRF-Chem 

• Compare the physical features of the simulated storm against 
aircraft and ground-based observations 

• Add flash rate parameterization schemes (FRPSs) to the model and 
identify the best match to observations 

• Determine NO production scenario for 1C and CG flashes following a 
lightning-generated NO x (LNO x ) scheme 
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Background 


• Storm system developed ~21Z May 29 
along KS/OK border and continued until 04Z 
May 30 

• Aircraft sampled storm and its environment 
from 20Z May 29 to 01Z May 30 

- DC-8 focused on storm inflow and outflow 

- GV and Falcon concentrated on outflow 


N EXRAD Composite Reflectivity 2240Z on 29 May 
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• Ground-based instrumentation included: 

— Dual-Doppler radar (NEXRAD level II 
regional; Data courtesy of C. Homeyer) 

- National Lightning Detection Network 
(NLDN) cloud-to-ground flash data 

- Oklahoma Lightning Mapping Array (LMA) 
flash initiation density data 



Blue circles: LMA stations 

Green outline: Extent of 3-D lightning mapping capability 
Gray outline: Extent of 2-D lightning detection 
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WRF-Chem Model V3.5 


• Nested domains: 15-km and 3-km 

• Initialized with DART and GFS for 
boundary conditions 

• Used coarsely prescribed IC:CG 
ratios (Boccippio et ai, 2001) 


LOCATION OF WRF PROCESSORS 
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+ Grid Cell Center LONGITUDE (°W) □ Processor 


Type of Scheme 

Selection for Simulation 

Microphysics 

Morrison 

Planetary boundary layer 

Yonsei University (YSU) 

Radiation 

Rapid radiative transfer model for GCMs (RRTMG) 

Flash rate 

Maximum vertical velocity (W max ) 

Lightning-generated NO x (LNO x ) 

DeCaria et al. (2000, 2005) 
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Flash Rate Parameterization Schemes 


• Based on simulated thunderstorm's physical features 

• Six types have previously been used in cloud- 
resolving models: 


Type of FRPS 

Equation (flashes min' 1 ) 

References 

Maximum vertical velocity 

5.7 x 10' 6 x W 4 5 

max 

Price & Rind, 1992 

Cloud top height 

3.44 x lO 5 x H 4 - 9 

Price & Rind, 1992 

Updraft volume 

6.75 x 10 11 x w 5 - 13.9 

Deierling & Petersen, 2008 

Ice water path 

33.33 x IWP-0.17 

Petersen et al., 2005 

Ice mass flux product 

9.0 x 10- 15 x (f p x f np ) + 13.4 

Deierling, 2006; 
Deierling et al., 2008 

Precipitation ice mass 

3.4 x IQ’ 8 x p m - 18.1 

Deierling et al., 2008 
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LNO x Parameterization Scheme 
(DeCaria et al v 2005) 

• Gaussian vertical distribution 
of 1C (bimodal) and CG (single 
mode) NO production based 
on typical lightning flash 
channel distributions 

• Lightning channels set to 
maximize at -15°C (CG and 1C) 
and -45°C (1C) 

• 500 moles NO per 1C and CG 
flash (Ott et a!., 2010) 

• Horizontal placement of NO 
based on reflectivity > 20 dBZ 
in each grid cell 


(o) 


CG flashes 


(b) 



1C flashes 



f(2) (km- 1 ) 


m (km-’) 
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Methodology 


Created moving spatial masks at 
10-min intervals for comparison of 
observed and model-simulated storms 

Used offline calculations, with 
adjustment factors, to analyze the 
six FRPS trends 

Calculated NLDN total flashes given 
NLDN CG flashes and mean IC:CG ratio 
for the storm region (3.9 ± 0.49), which 
is based on Boccippio et al. ( 2001 ) 



GOES-13 1-km Visible at 0008Z (NCAR/EOL) 


Total flashes = CG flashes x [1/0.93] NLDN DE x [IC:CG ratio + 1] 


Compared flash rate trends over the observed and model-simulated 
storm's lifetime 
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Initial Comparison of Storm Features 


Model-Simulated Composite Reflectivity 2300 Z on 30 May 
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Black rectangle represents the spatial mask surrounding the 
model-simulated cell of interest at 2300Z May 30 


NEXRAD Composite Reflectivity 0010Z on 30 May 
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Reflectivity (dBZ) 

Red rectangle represents the spatial mask surrounding the cell of 
interest at 0010Z May 30 

Model-simulated storm: 

• Began ~1-1.5 hour before observed storm 

• Exceeds area of observed storm by roughly a 
factor of two 
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Compared instantaneous 
flash rates from WRF at 
10-min intervals with 
corresponding 1-min 
periods from the observed 
NLDN flash rates 


Model-simulated flash 
rate trends are adjusted 
90 minutes later to 
coincide with observations 


Adj. instant 10-min WRF output vs. instant 10-min Adj. NLDN Obs. 
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Offline Calculation 
0 min = 2000Z 


■ Observations 

WRF - precipitation ice mass 
• WRF - updratt volume 
WRF ■ maximum vertical velocity 
WRF - ice mass flux product 
WRF - ice water path 
WRF - doud top height 


0 80 1 89 1B0 240 300 300 

Time (minutes) 


Flash Rate 

Parameterization Scheme 

Total Flashes 
Prior to Scaling 

Scaling Factor 

Maximum vertical velocity 

3,951 

1.1310 

Cloud top height 

708 

6.3138 

Updraft volume 

21,118 

0.2116 

Ice water path 

4,452 

1.0035 

Ice mass flux product 

36,745,336 

0.0001 

Precipitation ice mass 

164,749 

0.0271 


• Flux product, precipitation ice mass, and 
updraft volume trends are similar to the 
increasing trend of observations 

• Timing of W max and ice water path peaks is 
similar to observations (140, 200, & 310 min) 

• Magnitude of observed primary peak greater 
than those in FRPSs 

• W max and ice water path schemes need the 
least adjustment to match observed total 
flashes at each 10-min interval (4,468 
flashes) 
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* Plots courtesy of M. Bela 


Observations taken in cloud-free air to the 
south of the storm system 

To define the real background air at 9-12 
km, the 10 th percentile is used to remove 
any influence from old convective outflow 


Aircraft 

Time of UT Background 
Sampling (Z) 

DC-8 

20:40-21:10 (black) 

22:15-22:30 (red); 
22:58-23:05 (yellow) 
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Model-simulated vertical cross-section 
taken in cloud-free air to the south of the 
storm system 


Aircraft 

Time of UT Background 
Sampling (Z) 

Model 
Time (Z) 

DC-8 

20:40-21:10 

20:00 

22:15-22:30; 

22:58-23:05 

21:00; 

22:00 


Altitude 

NO 

(PPb) 

no 2 

(PPb) 

NO x 

(PPb) 

6-9 km 

Obs. 

0.025 

0.034 

0.059 

{Median) 

Model 

0.020 

0.017 

0.01-0.05 

9-12 km 

Obs. 

0.029 

0.033 

0.062 

( 10 th percentile) 

Model 

0.041 

0.023 

0.05-0.5 


Model-simulated NO x at 9-12 km is given 
as a range. Although the observed NO x is 
found at the lower end of this range, it 


should be kept in mind that the observed 
value represents the 10 th percentile. 
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Location of vertical cross-sections 



LNO x forecast for Oklahoma 
convection the morning of 
flight. Forecast based on WRF 
ARW 00Z May 29 model run. 

Will compare LNO x prediction 
of 2-3 ppbv in forecasted storm 
against observed aircraft 
measurements and WRF-Chem 
model-simulated NO x 


20 1 2 - 05-29 20 : 00 :00 
Lightning NOx Tracer (ppbv) 
OK-TX cross-sections 
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Initial LNO x Analysis from Aircraft 


N EXRAD Composite Reflectivity 2230Z on 29 May 



N EXRAD Composite Reflectivity 0040Z on 30 May 




Radar observations overlaid with 10-min 
intervals of aircraft NO x measurements 

NO x measurements increase from 
sampling in cloud-free upper tropospheric 
air (0.02-0.1 ppbv) to making transects 
through anvil outflow (peaks to ~2 ppbv) 


14 


Two methods for calculating NO x production per flash 


Molecules NO x estimated from volume 
Ridley et al. (1996, 2004) 


Top View 


*Volume= sum of surface areas from 
NEXRAD CAPPI images (2 km resolution) 



P(NO x ) oc LNO x enh * Volume 


P(NOJ in units of molecules 


NO x flux out of anvil 
Chameides et al. (1987), 
Huntrieser et al. (1998, 2002) 

Vertical Cross Section 
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anvil depth 
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P(NOJ in units of molecules s 
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Divide by #flashes s -1 to get molecules flash -1 


Divide by #flashes to get molecules flash -1 
*Slide courtesy of I . Pollock 



Total number of lightning flashes during the 29-30 May 2012 storm: 


Storm start 

Storm end 

Duration 

CG flashes 

Total flashes (NLDN) 

Total flashes 

(UTC) 

(UTC) 

(hrs) 

(NLDN) 

IC:CG=2 

IC:CG=10 

(LMA) 

29 May 
21:00 

30 May 
00:49 

3.82 

2,851 

8,553 

64,361 

69,650 


Using volume method to find P(NO x ) in units of (xlO 25 molecules flash _1 ): 


Mean LNOx 
(xlO 15 molec nrr 3 ) 

Storm volume 
(xlO 13 m 3 ) 

P(NO x )nldn 

P(NO x )lma 

30 dBZ 20 dBZ lOdBZ 

30 dBZ 

20 dBZ 10 dBZ 

30 dBZ 

20 dBZ 10 dBZ 



3.7 


0.5 


5.2 ±0.1 

6.0 16.2 40.6 

61.4 

(mol flash 1 ) 

9.9 24.8 

8.3 

(mol flash 1 ) 

1.2 3.0 


Using flux method to find P(NO x ) in units of (xlO 25 molecules flash -1 ): 


Mean LNO x 
(xlO 25 molec s -1 ) 

NLDN flash rate 
(flashes s -1 ) 

P(MO x )nldn 

LMA flash rate 
(flashes s -1 ) 

P(NO x )lma 



4.2 


0.5 

2.7 ±0.1 

0.64 

69.7 

5.20 

8.3 



(mol flash 1 ) 


(mol flash 1 ) 


*Slide courtesy of I. Pollock 

















Conclusions 


• Based on offline calculations, W max FRPS was selected for use in model: 

- Needs little adjustment to match the observed total flashes 

- Coincides with several of the observed flash rate peaks 

• Scale up model-simulated flash rates in offline calculations and scale down online: 

- May partly be due to how offline calculations are computed 

• Model overestimate of observed flashes may be due to: 

- Area of model-simulated storm ~2x larger than observed 

- Observed storm passes over northern edge of LMA 

• Initial look at NO x chemistry in UT air undisturbed by storm: 

- At 6-9 km, NO values similar between aircraft and model, and model-simulated N02 
underestimates observations by ~0.02 ppbv 

— At 9-12 km, the 10 th percentile NO x values are similar (~0.06 ppbv) between the aircraft and 
model 

• Generate similar P(NO x ) when using volume (> 30 d BZ) and flux methods 

- Estimated P(NO x ) is much smaller than 500 moles flash 1 used in WRF-Chem 

• WRF-Chem model estimates of P(NO x ) in works 
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Future Work 


• Perform a trace gas simulation 
and analysis of NO x , CO, and 
0 3 using WRF-Chem 

• Compare model-simulated 
LNO x against aircraft 
measured NO x 



• Determine NO production 
scenario per 1C and CG flash 
that best matches aircraft 
observed NO x mixing ratios 



• Investigate 0 3 changes 
downwind of flight 
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